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The physically interesting gravitational analogue of magnetic monopole in electrodynamics is 
considered in the present paper. The author investigates the field equation of gravitomagnetic 
matter, and the exact static cylindrically symmetric solution of field equation as well as the motion 
of gravitomagnetic charge in gravitational fields. Use is made of the mechanism of gravitational 
Meissner effect, a potential interpretation of anomalous, constant, acceleration acting on the Pioneer 
10/11, Galileo and Ulysses spacecrafts is also suggested. 
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I. INTRODUCTION 

It is well known that the field equation of general relativity under low-motion weak-field approximation is anal- 
ogous to Maxwell's equation of electromagnetic field. This similarity leads us to consider the following interesting 
gravitational analogues of electromagnetic phenomena: (i) in electrodynamics a charged particle is acted upon by 
Lorentz magnetic force, in the same manner, a particle is also acted upon by the gravitational Lorentz force in 
weak-gravity theory [1,2]. According to the principle of equivalence, further analysis shows that in the non-inertial 
rotating reference frame, this gravitational Lorentz force is just the fictitious Coriolis force [3]; (ii) there exists the 
Aharonov-Bohm effect in electrodynamics [4], accordingly, the so-called gravitational Aharonov-Bohm effect, i.e., the 
gravitational analogue of Aharonov-Bohm effect also exists in theory of gravitation, which is now termed Aharonov- 
Carmi effect [5-7] ; (iii) a particle with intrinsic spin possesses a gravitomagnetic moment of such magnitude that it 
equals the spin of this particle which can be coupled to gravitomagnetic fields. Mashhooms spin-rotation coupling is 
the non-inertial realization of the interaction of gravitomagnetic moment with gravitomagnetic field [8,9]. We think 
that investigation of general characters of gauge field deserves generalization to that of the gravitational field [10,11]. 
In this paper, we study another physically interesting phenomenon, i.e., the gravitational analogue of magnetic charge 
in electrodynamics. 

In electrodynamics, electric charge is a Noether charge while its dual charge ( magnetic charge ) is a topological 
charge, since the latter is related to the singularity of non-analytical magnetic vector potentials. Magnetic monopole 
[12] attracts attentions of many physicists in various fields such as gauge field theory, grand unified theory, particle 
physics and cosmology [13-17]. In the similar fashion, it is also interesting to consider gravitomagnetic charge which 
is the source of gravitomagnetic field just as mass ( gravitoelectric charge ) is the source of gravitoelectric field ( 
Newtonian gravitational field ). In this sense, gravitomagnetic charge is also termed dual mass. It should be noted 
that the concept of mass is of no significance for the gravitomagnetic charge, and then it is of interest to investigate the 
relativistic dynamics and gravitational effects of this topological dual mass. Historically, gravitomagnetic monopole 
attracts attentions of many authors in both classical dynamics and quantization formulation [18-23]. For historical 
review one may be referred to the paper by Linden-Bell and Nouri-Zonoz [24]. 

In literature, to the best of our knowledge, many authors were concerned with monopole rather than with charge, 
namely, they considered equations of motion of gravitomagnetic monopoles more than field equations of gravitomag- 
netic charges/matter. In this paper, the author proposes the field equation of gravitomagnetic matter and investigates 
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the equation of motion of gravitomagnctic monopole, and then exactly solves the static cylindrically symmetric so- 
lution and discusses the motion of a photon in the gravitomagnctic fields. In Sec. V, a potential interpretation of 
anomalous, constant, acceleration acting on the Pioneer 10/11, Galileo and Ulysses spacecraft [25] is suggested. The 
relation between the non-analytic property of metric and the topological gravitomagnetic charge is discussed in Sec. 
VI. In Sec. VII, the author concludes with some remarks. 



II. TOPOLOGICAL CHARGE AND FIELD EQUATION OF GRAVITOMAGNETIC MATTER 

First we take into consideration the 5*0(3) non-Abel magnetic monopole proposed by Li et al. [26,27], which enables 
one to consider the topological dual mass in curved spacetime, then we suggest the gravitational field equation of 
gravitomagnetic charge. 

The linear element on the two-dimensional spherical surface is ds 2 = r 2 d9 2 + r 2 sin 2 9dtp 2 where the non-vanishing 
Christoffel symbols T % ^- (y(i, i, j = 1, 2) is given as follows 

r^ = r 2 2 = ^, ri 2 = -sintfcostf (2.1) 

and its matrix form (Tg)j and (r v )* are therefore 

fO \ / -sinflcosfl\ 

Ue) = n cose , (M = cose n • (2-2) 

V u sine / V sine u / 

Under the transformation M" 1 , the basic vector rsm9dcp that varies with 9 may be transformed into what is 9- 
independent, namely, 

dip J \r sin 9d(p J \ r sin 9 J 

Accordingly, is transformed into = MT^M -1 + Mc^M -1 and the results are of the form 

-cc 
v cos 9 

with Xs = ^ ^ q ^ being one of the generators of 50(3) group. There exists the connection between f M and the 

potential which may be written as = gW^X 3 and we therefore obtain gWg — 0, gW^ = — cos 9, and the 
non- vanishing field strengths are written in the forms 

,3 _ ,3 _ d K dW$ _sm6 

Je v ~ Uo- d6 d V ~ g ■ [ ' 

It follows that the source charge g of field strength is expressed by 



r = o, t v = [ 2, a T d )=- cosdx z ( 2 - 4 ) 



fl v d6dip=j. (2.6) 



It should be noted that g is Noether charge while g is topological charge whose existence is related to the singularities 
of field potentials. The work of Li et al. [26] leads us to consider the dual charge associated with spacetime. In fact, 
although no topological charge g exists seen from in the three-dimensional space, the reason why the two-dimensional 
surface is curved seen from the observer on the above two-dimensional spherical surface can be equivalcntly ascribed 
to the presence of this kind of dual charge. The concept of above topological charge may be readily generalized to 
that of gravitomagnetic charge. 

The dual of Riemann curvature tensor, R^js, may be defined as follows 

RpvyS — tfn/*^ R^Saf) + tys"^ Rfivafl (2-7) 

with e A1!y Q/3 being the completely antisymmetric Levi-Civita tensor satisfying its covariant derivative e^ u a = 0. 
It should be noted that R^jS = in the absence of gravitomagnetic matter since no singularities associated with 
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topological charge exist in the metric functions and therefore Ricci identity still holds. However, once the metric 
functions possess non-analytic properties in the presence of gravitomagnetic matter ( should such exist ), the dual 
curvature tensor is therefore no longer vanishing due to the violation of Ricci identity. The dual curvature scalar that 
contains only the first derivative and second derivative of field g^ v can be chosen to be R = — e^" 13 R^a/s ■ By making 
use of the variational principle, the dual Einstein's tensor is obtained as follows 

s I -V=9^ af3 R^apdn = \ [ V^g(R„, - ^ „ a(3 R„ Safi )6g^<m 

Jn 1 Jn 

= I V^gG^Sg^dn (2.8) 
Jn 

with R^ = g' JT R <J fiuT- Further calculation yields 

Gfi V = ^ R-naf3-y ~ e ^ Rvapf), (2-9) 

which is considered the dual of the Einstein's tensor arising on the left handed side of gravitational field equation of 
gravitomagnetic matter. Note, however, that once the gravitomagnetic charge is absent in spacetime, vanishes 
due to the Ricci identity. But, the non- analytic metric functions caused by the presence of gravitomagnetic charge 
may result in the violation of Ricci identity and therefore lead to G^ v 7^ 0. Since the dual Einstein's tensor is an 
antisymmetric tensor, in what follows we will construct the antisymmetric source tensor of gravitomagnetic charge 
that is on the right handed side of gravitational field equation. We show that for the Fermion field, the antisymmetric 
tensors constructed in terms of field tp and space-time derivatives are of the form 

K^ = i^d v - 7v d^, H liV = € ltu a ^K aP (2.10) 

and regard the linear combination of them as the source tensor in the field equation of gravitomagnetic charge, where 
7^s denote the general Dirac matrices with respect to and satisfy 7^7„ + 7^7^ = ^g^v Thus the field equation 
governing the distribution of gravitation of gravitomagnetic charge may be given as follows 

= KiK^ + k 2 H^ (2.11) 

with K\,K2 being the coupling coefficients between gravitomagnetic matter and gravity. One of the advantages of 
this field equation is that it does not introduce extra tensor potentials when allowing for gravitomagnetic monopolc 
densities and currents. This fact is in analogy with that in electrodynamics, where the equation d^F^" = governs 
the motion of electromagnetic fields produced by magnetic monopoles. Additionally, further analysis shows that the 
linear combination of the following antisymmetric tensors, — JvJn)^ an d i^^^i'ilalf} — Ipla)^ serves as 

the cosmological term in the above gravitational field equation of Fermion field. 

It is believed that there would exist formation ( and creation ) mechanism of gravitomagnetic charge in the grav- 
itational interaction, just as some prevalent theories [16] provide the theoretical mechanism of existence of magnetic 
monopole in various gauge interactions. Magnetic monopole in electrodynamics and gauge field theory has been 
discussed and sought after for decades, and the existence of the 't Hooft-Polyakov monopole solution has spurred new 
interest of both theorists and experimentalists [16,28,29]. Similar to magnetic monopole, gravitomagnetic charge is 
believed to give rise to such situations. If it is indeed present in universe, it will also lead to significant consequences in 
astrophysics and cosmology. We emphasize that although it is the classical solution to the field equation as discussed 
above, this kind of topological gravitomagnetic monopole may arise not as fundamental entities in gravity theory. 



III. MOTION OF GRAVITOMAGNETIC CHARGE 



According to its gravitational properties, gravitomagnetic charge can be called dual mass, and accordingly, a point- 
like particle possessing mass is also called gravitoelectric monopole. This, therefore, implies that the concept of mass 
is of no essential significance for gravitomagnetic matter; it is of interest to investigate the motion of gravitomagnetic 
monopole in curved spacetime. Although Ricci identity is violated due to the non-analytic properties caused by the 
existence of gravitomagnetic charge, Bianchi identity still holds in the presence of gravitomagnetic charge. It follows 
that the covariant divergence of G^ vanishes, namely, 

G^. u = 0. (3.1) 
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Then in terms of the following field equation 



= S" v (3.2) 

with the antisymmetric source tensor of gravitomagnetic matter S^ u being KyK^ + n 2 H^ v , one can arrive at 

S^. v = (3.3) 

which can be regarded as the equation of motion of gravitomagnetic charge in the curved spacetime. It is useful to 
obtain the low-motion and weak- field-approximation form of Eq. (3.3), which enables us to guarantee that Eq. (3.3) 
is indeed the equation of motion of gravitomagnetic monopole. 

The general Dirac matrices under the weak- field approximation may be obtained via the relations ^^v+lvl^ = ^g^u 
and the results are given as follows 



(3.4) 



where i = 1,2,3; /3 = 7^. 7^ and 7^ are the constant Dirac matrices in the flat Minkowski spacetime. The 



— , and gravitomagnetic vector potentials are g % = g 



gravitoelectric potential is defined to be g° = 
In the framework of dynamics of point-like particle, the source tensor is therefore rewritten as 

S"" = p [k^XF - g v U») + K 2 e^(g a U fj - g fj U a )} , 



Q-i 



1,2,3). 



(3.5) 



where p denotes the density of gravitomagnetic matter. It follows from Eq. (3.3) and Eq. (3.5) that there exists the 
gravitational Lorentz force density in the expression for the force acting on the gravitomagnetic charge, namely, by 
ignoring some small term and using the static condition gfnS = 0, one can obtain the following expression ( in the 
unit c = 1 ) 



)]- 2K2gX (^ V + V ' 9 °) 



+Ki0°^og-2 K25 ° (V. 9 ° x v) 



(3.6) 

with v being the velocity of the tested gravitomagnetic monopole. It is apparent that V x g — v x (V.g° — gfrrg) is the 
expression associated with gravitational Lorentz force density. Note that in Eq. (3.6) k\,k 2 are considered coupling 
constants. However, further analysis shows that at least one of them is not a constant and if the relation 



= 2k 2 

between them is assumed, then Eq. (3.6) may be rewritten as 

_d_ 



V x g - v x ( V fl ° - ^g 



(3.7) 



(3.8) 



where we ignore the small term of second order and the derivative term of coupling coefficients with respect to 
spacetime coordinate, . It is well known that the form of Eq. (3.8) is the equation of motion of a particle 
acted upon by the Lorentz force. Hence, Eq. (3.3) is believed to be the generally relativistic equation of motion of 
gravitomagnetic monopole in the Ricmann spacetime. 

Investigation of relativistic dynamics of the topological dual mass is of interest. It should be noted that although 
gravitomagnetic monopole does not possess mass, it still has energy. Since the dual mass is a kind of topological charge 
which is very different from Noether charge, no mechanism of interaction may turn it into mass (gravitoelectric charge), 
and vice versa. 

The gravitational analogue of Meissner effect in superconductivity is gravitational Meissner effect. Due to the 
conservation law of momentum, mass-current density may be conserved in the process of scattering in perfect fluid, 
which is analogous to the superconductivity of superconducting electrons in superconductors cooled below T c . Since 
gravitational field equation under linear approximation is similar to the London s equations of superconductivity, one 
can predict that gravitational Meissner effect arises in perfect fluid. The author holds that the investigation of both 
the effect of gravitomagnetic matter and gravitational Meissner effect may provide us with a valuable insight into the 
problem of cosmological constant and vacuum gravity [30-33]: the gravitoelectric field ( Newtonian field of gravity 
) produced by the gravitoelectric charge ( mass ) of the vacuum quantum fluctuations is exactly cancelled by the 
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gravitoelectric field due to the induced current of the gravitomagnetic charge of the vacuum quantum fluctuations; the 
gravitomagnetic field produced by the gravitomagnetic charge ( dual mass ) of the vacuum quantum fluctuations is 
exactly cancelled by the gravitomagnetic field due to the induced current of the gravitoelectric charge ( mass current 
) of the vacuum quantum fluctuations. Thus, at least in the framework of weak-field approximation, the extreme 
space-time curvature of vacuum caused by its large energy density does not arise, and the gravitational effects of 
cosmological constant is eliminated by the contributions of the gravitomagnetic charge ( dual mass ) . If gravitational 
Meissncr effect is of really physical significance, then it is necessary to apply this effect to the early universe. Some 
related topics such as gravitational Hall effect and gravitational magnctohydrodynamics may be considered for further 
consideration. 



IV. STATIC CYLINDRICALLY SYMMETRIC EXACT SOLUTION OF FIELD EQUATION 

In this section the static cylindrically symmetric gravitomagnetic field and the evolution of wavefunction of photon 
in gravitomagnetic field are considered. Suppose that the form of linear element describing the static cylindrically 
symmetric gravitomagnetic field is given by 

ds 2 = d(a; ) 2 - Ax 2 - dy 2 - dz 2 + 2g 0x {y)dx°dx + 2g 0y {x)dx°dy, (4.1) 

where we assume that the gravitomagnetic potentials go x and goy are functions with respect to y and x, respectively. 
Thus we obtain all the only nonvanishing values of Christoffel symbols as follows: 



r -r - 1 ( d90x i d 9Qy \ r _ r _ _ r _ _ r _\( d 9ox dg 0y \ 

L °' xv ~ L °' yx ~ 2\ d y + dx J ' L *>°y- L *'V°- L v> *- L y' x0 ~2\dy 8x ) 



(4.2) 



Since the field equation of gravitomagnetic matter is the antisymmetric equation, we might as well take into account 
a simple case of the following simplified equation 

e 0a ^R° aPl = P m (4.3) 

with pm being the parameter associated with the coupling constants and gravitomagnetic charge. It is therefore 
apparent that Eq. (4.3) agrees with Eq. (3.2). Hence, the solution of the former equation also satisfies the latter. 
For the reason of the completely antisymmetric property of the Levi-Civita tensor, the contravariant indices a, (3, 7 
should be respectively taken to be x,y, z, namely, we have 

e 0a ^R° aPl = 2e te »*(i2° xyz + R° zxy + R° yzx ). (4.4) 

There exist the products of two Christoffel symbols, i.e., g aT (T T , Q7 I\ CT( 3 — T T . a pT\, ai ) in the definition of the Riemann 
curvature, R\ a p-y- Apparently, the products of two Christoffel symbols ( i.e., the nonlinear terms of field equation ) 
contain the total indices, x, y, z of three-dimensional space coordinate ( namely, these indices are taken the permu- 
tations of x,y,z ) and therefore vanish, in the light of the fact that the Christoffel symbol with index z is vanishing 
according to Eq. (4.2). 

In view of the above discussion, one can conclude that Eq. (4.3) can be exactly reduced to a linear equation. 
It is easily verified that R\ a (3-y (A = x, y, z) vanishes with the help of the linear expression for R\ a f}j given by 

Rxafi-t = \ ( dJalfl + dx^dxi - dx^dx-y ~ dx^dxU ) and thc linear element expressed by Eq. (4.1). We thus obtain 
that J?° af3j — 3 00 i?oa/37- By the aid of thc following expression 

one can arrive at 

r=9 



e^RO a/37 = -4=V-(Vxg), (4.6) 



where the gravitomagnetic vector potentials, g, are defined to be g = (—gox,—goy,—goz)- Substitution of Eq. (4.6) 
into Eq. (4.3) yields 
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V • (V x g) 



-g 

-WP M - 



(4.7) 



Note that Eq. (4.7) is the exact static gravitational field equation of gravitomagnetic matter derived from Eq. (3.2), 
where use is made of the expression (4.1) of linear element. For the case of cylindrically symmetric gravitomagnetic 
field with the nonvanishing pu being present only in the x-y plane with z — 0, the gravitomagnetic field, B g , which 
is defined as V x g, may be written 



(B fl ), = - 



cm 



-9 



,,00 



(4.8) 



2 = 



with (7m being the surface density associated with pu and ( ] denoting the value of in the x-y plane 

V 9 J 2=0 9 

where z = 0. It follows from Eq. (4.8) that the direction of B g is parallel to the z-axis. The metric components, 
9ox,9oy, are therefore readily obtained as follows 



90x = -yt/, 



90y = — 2"&> 



(4.9) 



With Bg = -"-f- „ 

In order to obtain the contravariant metric g^ v , we calculate the inverse matrix of the metric {g^ v ) and the result 
is given as follows 



or) 



4 



/ 1 




5oy 





\ 


9ox 


-(l+9o 2 y ) 


90x90y 







90y 


goxgoy 









V o 





-(14 


5te 4 


4/)/ 



(4.10) 



It is well known that a photon propagating inside the noncoplanarly curved optical fiber that is wound smoothly 
on a large enough diameter [34-36] is acted upon by an effective Lorentz force which may be expressed as (in the unit 
h = c = 1) [37] 



f = k= k x 



k xk 

k 2 



(4.11) 



where the effective magnetic field is with dot denoting the time rate of change of k(i). Eq. (3.8) has been 
shown to be the expression for the gravitational force acting on the gravitomagnetic monopole, in the similar manner, 
one can consider the motion of gravitoelectric charge, for instance, a photon propagating in the static cylindrically 
symmetric gravitomagnetic field. For the sake of analyzing this problem conveniently, we first take into account the 
time evolution of wavefunction of a photon in the weak gravitational field. The infinitesimal rotation operator of 
wavefunction of the photon in gravitomagnetic field is given by Ur — 1 — «Ai9 • J with [26] 



k(t) x [k(t) + Ak] _ k(i)x k (t) 



fc 2 



fc 2 



At, 



(4.12) 



where Ak is defined to be k At. Simple calculation shows that the effective Hamiltonian is of the form H e ff = - J. 
Given that this infinitesimal rotation of wavefunction of the photon is caused by the gravitomagnetic field, it follows 
from both Eq. (4.11) and (4.12) that the equation of motion of a photon in gravitomagnetic field may be written in 
the following form 



k xk 

k 2 



= V x g. 



(4.13) 



This formulation is readily generalized to the case of massive particle moving in gravitomagnetic field and a number 
of related topics concerning the matter wave in gravitomagnetic field may be further investigated. 
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V. A POTENTIAL INTERPRETATION OF ANOMALOUS ATTRACTIVE FORCE ACTING ON 

PIONEER SPACECRAFTS 



Taking the effects of gravitomagnetic charge into consideration is believed to be of essential significance in resolving 
some problems and paradoxes. An illustrative example that has been briefly discussed is its application to the 
problem of cosmological constant. Additionally, in 1998, Anderson et al. reported that, by ruling out a number 
of nongravitational potential causes such as the solar radiation pressure, precessional attitude-control maneuvers, 
nonisotropic thermal radiation, radiation of the spacecraft radio beam and so on, radio metric data from the Pioneer 
10/11, Galileo and Ulysses spacecrafts indicate an apparent anomalous, constant, acceleration acting on the spacecraft 
with a magnitude ~ 8.5 x 10~ 8 cm/s 2 directed towards the Sun [25]. Is it the effects of dark matter or a modification 
of gravity? Unfortunately, neither easily works. It is interesting that, by taking the cosmic mass, M = 10 53 kg, and 
cosmic length scale, R = 10 26 m, our calculation shows that this anomalous acceleration is just equal to the value of 
gravitational field strength on the cosmic boundary due to the total cosmic mass. This fact leads us to consider a 
theoretical mechanism to interpret this anomalous phenomenon. The author favors that the gravitational Meissncr 
effect may serve as a potential interpretation. Here we give a rough analysis, which contains only the most important 
features rather than the precise details of this theoretical explanation. Once gravitomagnetic matter exists in the 
universe, parallel to London's electrodynamics of superconductivity, gravitational field may give rise to an effective 
mass Trig = -^yJ&TrGpm due to the self-induced charge current [38], where p m is the mass density of the universe (For 
the case of arbitrarily strong gravitational fields, see the interesting work of Visser [39] , where he defined the mass 
term of graviton by introducing a non-dynamical background metric). A constant acceleration, a, may result from 
the Yukawa potential and can be written as 



a = 



GM (m a c\ 2 GM 

MPmG) 
'■M 

(5.1) 



GM G(4irR 3 p m ) GM 



c 2 R R 2 ' R 2 ' 

where for the universe, use is made of ~ 0(1), ^R? pm — M, which holds when approximate estimation is 
performed. Note, however, that this is an acceleration of repulsive force directed, roughly speaking, from the center 
of the universe. By analyzing the NASA'S Viking ranging data, Anderson, Laing, Lau et al. concluded that the 
anomalous acceleration does not act on the body of large mass such as the Earth and Mars. If gravitational Meissner 
effect only affected the gravitating body of large mass or large scale rather than spacecraft, then seen from the Sun 
or Earth, there exists an added attractive force acting on the spacecraft. In the following, gf roc denotes the free-body 
acceleration caused by the cosmic mass, and apparently the velocity of the Sun and the spacecraft agrees with the 
following equation (using approximate analysis) 

d_ d_ 

~j^" v sun Sfrcc T" 3-; v spacccraft — Sfrcci (5-2) 

then the velocity of the spacecraft relative to the Sun satisfies 

(^spacecraft — v sun) = a , (35) 

where the acceleration of spacecraft due to the Solar gravitation has been ignored. The negative sign on the right 
handed side of Eq. (35) shows that this added force gives rise to an anomalous, constant, acceleration directed towards 
the Sun. 

Anderson et al. reported that without using the apparent acceleration, i.e., using an independent analysis, the 
compact high accuracy satellite motion program (CHASMP) shows a steady frequency drift of about —6 x 10~ 9 Hz/s, 
or 1.5 Hz over 8 year. This equates to a clock acceleration, at of —2.8 x 1CU 18 s 2 /s. If there were a systematic drift 
in the atomic clocks of the DSN (Deep Space Network) or in the time-reference standard signals, this would appear 
like a nonuniformity of time; i.e., all clocks would be changing with a constant acceleration. Anderson et al. believed 
that people cannot rule out this possibility without adequate analysis. However, we think that the so-called clock 
acceleration has no fundamental significance, but favors the existence of the anomalous acceleration acting on these 
spacecraft, seen from the Sun. The relation between coordinate time dx° in the solar system and the proper time dr 
is 

dx° (j) sun + ar 

^- l -^ 2 ~ ^ 
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with sun being the gravitational potential caused by the mass of the Sun and r the radial distance from the Sun. By 
ignoring ^p 2 - that is familiar to us, the added clock acceleration may be written as 

d 2 x° a dr a 
°* cdr 2 c 2 dr c ' 

where ^ is taken to be c, the speed of the radio signals from the spacecrafts. It follows from Eq. (5.4) that 

a = —cat- (5-5) 

This expression can be proved correct by inserting the experimental values, a t — —2.8 x 10~ 18 s 2 /s and a = 8.5x 
10~ 10 m/s 2 , into the above equation. 

Although it seems natural to interpret this anomalous gravitational phenomenon, there is still something that 
deserves further research. Why does not the gravitational Meissner effect explicitly affect the body of small mass? 
Maybe small-mass flow cannot serve as the self-induced charge current that provides the gravitational field with 
an effective mass in the field equation. The theoretical resolution of this problem is not very definite at present. 
However, the sole reason that the above resolution of the anomalous acceleration is somewhat satisfactory lies in 
that no adjustable parameters exist in this theoretical mechanism. It is one of the most important advantages in the 
above mechanism of gravitational Meissner effect, compared with some possible theories of modification of gravity 
[40], which are always involving several parameters that cannot be determined by theory itself. These theories of 
modification of gravity were applied to the problem of the anomalous acceleration but could not calculate the value of 
the anomalous acceleration consistent with experiments performed by Anderson et al. . In spite of some defects, the 
theoretical mechanism of gravitational Meissner effect suggested above still gives a valuable insight into the problem 
of anomalous acceleration acting on the Pioneer spacecrafts. 



VI. NON-ANALYTIC PROPERTY OF METRIC AND TOPOLOGICAL GRAVITOMAGNETIC CHARGE 

Since the gravitomagnetic charge is topological charge of spacetime, it may lead to the non-analytic property of 
space-time metric, and result in the violation of Ricci identity. This may be illustrated as follows: 

Based on the gravitational field equation (2.11) of gravitomagnetic charge, we calculate the following formula 

e af3 ^Gfxv = e ap ^{mK^ + k 2 H^), (6.1) 

and obtain 

€ af3 ^"(KiK^ + K 2 H^) = K!H a p + 4K 2 K a/3 , 

z a p ""GW = 2{R al3 - R 0a ), (6.2) 

or 

Rap - Rf3a = -^KlH a fj + 2K 2 K aj 3, (6.3) 

where use is made of e a/ 3i^e AcrrAI = .9^5^.9^] with A, a, r and a, (5, v being completely antisymmetric, respectively. It 
follows from Eq. (2.11) and (6.3) that 

2n 2 {R a fj - Rpa) ~ KxG a {i = (4/t 2 - Kl)K a/ 3 7 

2Ki(R a p - Rp a ) - 4K 2 G a p = {k\ - 4 K 2 2 )H a/3 . (6.4) 

It is known that R a p = Rp a holds in the absence of gravitomagnetic charge in general relativity. However, when there 
exists gravitomagnetic charge, the symmetric property of Ricci tensor is no longer valid. Further calculation shows 
that 

R a0 -R 0a -^9 [q^^-q^q^)' (6-5) 

and then it is apparently seen that the non-analytic property of metric leads to the non-integral condition, and that 
the above case is in exact analogy with that of the Aharonov-Bohm effect in electrodynamics. This, therefore, implies 
that gravitomagnetic charge possesses the topological or global properties, and that the geometric and non-analytic 
properties of spacetime deserve detailed study. 
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VII. CONCLUDING REMARKS 



The present paper considers the concept of gravitomagnetic charge and its dynamics in curved spacetime, con- 
taining the field equation and the equation of motion and the static exact solution. Differing from the symmetric 
property (with respect to the space-time coordinate) of gravitational field equation of gravitoelectric matter, the field 
equation of gravitomagnetic matter possesses the antisymmetric property. This, therefore, implies that the number of 
the non- analytic metric functions is no more than 6. Although we have no observational evidences for the existence 
of gravitomagnetic charge, it is still of theoretical significance to investigate the gravity theory of the topological 
dual mass. It is interesting to speculate on the possibility that the formation mechanism of gravitomagnetic charge 
may be found in the theory of gravitational interaction, just as the case of the magnetic monopole in various gauge 
interactions. Once it is present in universe indeed, gravitomagnetic charge leads to physically interesting consequences 
in astrophysics and cosmology, particularly in the physics of early universe. In the gravitational fields where grav- 
itomagnetic charge exists, the influence of gravitational Meissner effect on the cosmological evolution should not be 
neglected. In the regime of strong gravity such as the stages of quantum cosmology and inflationary model, gravita- 
tional Meissner effect may cause important effects. From the point of view of the framework of classical field theory, 
matter may be classified into two categories: gravitomagnetic matter and gravitoelectric matter, according to their 
different gravitational features 1 . Although it is the classical solution to the field equation, this kind of topological 
gravitomagnetic monopoles may arise not as fundamental entities in gravity theory. 

With foreseeable improvements in detecting and measuring technology, it is possible for us to investigate quantum 
mechanics in weak-gravitational fields [41,42] associated with gravitomagnetic fields, and the investigation of dual 
mass is therefore essential in gravity theory. A potential application of gravitational Meissner effect to the problem 
of the anomalous acceleration is of interest. Although there exists the curiosity as to whether the gravitational 
Meissner effect is universal for all the gravitating body, regardless of its mass and scale, it gives the theoretical value 
of anomalous acceleration which is somewhat consistent with the experimental value. This theoretical possibility is of 
use, given that we have no plausible explanation so far. It is apparently worthwhile to improve this mechanism based 
on the gravitational Meissner effect. 
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